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ABSTRACT 



Space-based projects are providing a wealth of high-quality asteroseismic 
data, including frequencies for a large number of stars showing solar-like os- 
cillations. These data open the prospect for precise determinations of key stellar 
parameters, of particular value to the study of extra-solar planetary systems. 
Given the quantity of the available and expected data it is important to develop 
efficient and reliable techniques for analyzing them, including the determination 
of stellar parameters from the observed frequencies. Here we present the SEEK 
package developed for the analysis of asteroseismic data from the Kepler mission. 
A central goal of the package is to obtain a fast and automatic determination of 
the stellar radius and other parameters, in a form that is statistically well-defined. 
The algorithms are tested by comparing the results of the analysis with indepen- 
dent measurements of stellar radius and mass, for a sample of well-observed stars. 
We conclude that the SEEK package fixes stellar parameters with accuracy and 
precision. 

Subject headings: methods: statistical - stars: fundamental parameters - stars: 
oscillations - stars: solar-type 
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ASTROPHYSICAL CONTEXT 



Solar-like oscillation occurs in stars with convective envelopes and are seen as acoustic 
modes (p modes) with low angul ar degree / < 3 and i ntern iediate to high radia l order 



10 < n < 30. After early work by 



Ando fc Osakil fllQTSf ) and 



Goldreich fc KeelevI fll977al Jbh. 



evidence has been found that the excitation of solar-like oscillation is caus ed by turbulent 



convection in the u 


ppermost part o 


" the € 


Goode et al. 


1992 




Goldreich et al. 


1994) 



' the envelopes of stars (see for example 



Balmforth 



1992 



1994j ). An important aspect of the driving caused by 
turbulent convection is that the excitation of the mode occurs at random times, and hence 
the process is known as stochastic excitation. One of the effects of this stochastic behavior is 
that the eigenmode phases are changing with time. In the Sun, the mode s with the largest 



amplitudes, arounc 



Chaplin et al. 



20021 : 



3000 /iHz, have a coherence time of less than 12 days (iGelly et al.l 
19971 ). We can contrast this value to classical pulsators where the coherence 
time is of the order of the stellar evolution time scale. Short but recent introductions to the 
topic of stochas tic mode s excit ation , and references to m odern approaches to that problem. 



can be found in 



Houdekl (120061 ) and 



Samadi et al 



(120081). 



The first star other than the Sun showing evidence for individu al frequencies of 
solar-like oscillation, rj Bootis, was observed by 



Kjeldsen et al. 



( 1l995l ). using a technique 



measuring the equivalent widths of the stellar Balmer lines over time. Most subsequent 
observations of solar-like oscillation have been made with Doppler-velocity measurements 
using a laboratory reference (typically an iodine cell) to measure the shift in a stellar 
spectrum. This technique fi nally permitted the first unambiguous detection of solar-like 



oscillation in a Centauri A (IBouchy fc Carrier 



2001 



2OO2I ). Since then, p-mode excitation 



has been detected with ground-based observations in more than a do zen main-sequence 



and su bgia nt stars. Lists of the se solar-like oscillators can be found in 



(2008 


) and 


Bruntt et al. 


(2010) 



Bedding fc Kjeldsen 
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Doppler-velocity measurement is the technique of choice to observe solar-hke oscillation, 
providing the least interference from stellar background contributions to the 'noise' in 
the observations. The ratio between the signal and the stellar background noise, that is 
non-periodic signals coming, for instance, from granulation, is far higher in Doppler-velocity 
than in photometric intensity measurements. Solar-like oscillations can also be observed in 
photometric intensity measurement, albeit with a higher stellar background noise, but with 
the advantage of allowing simultaneous observation of a large number of stars. However, 
owing to noise induc ed by the Ea rth's atmosphere such observations are essentially only 



possible from space ([Harvey 



19881 ). Space-based photomet ric data for solar-li ke oscillations 



have started to be o btained with the help of the CoRoT (IBaglin et al. 



(IBorucki et al. 



Kepler satellite ( 



20061 ) and Kepler 



2010f) missions. Resu lts coming from these satellites, in particular the 



Gilliland et al. 



20101 ). are indeed a motivation for the work undertaken 
in the present study. The development of an objective technique to interpret results from 
asteroseismology, more than being intrinsically desirable, is now made necessary by the era 
of high-speed space photometry, where data bases containing large numbers of light curves 
ready for asteroseismic study are made available. By objective asteroseismology we mean 
a technique that is not directed at every step by subjective choice of a scientist and that 
results in an estimate of global stellar parameters along with realistic, if not completely 
objective, uncertainties on these parameters. The need for efficient and reliable, automatic 
analysis of as teroseismic data will be further emphasized with the ESA PLATO mission 



(ICatala 



20091 ) which, if finally selected, will be launched around 2018. 



The Sun is evidently the canoni cal example of solar-lik e oscil 



in the review on helioseismology by 



Christensen-Dalsgaardl (120021 ) the study of the Sun's 



ation. As shown in detail 



oscillations has proven that physical information can be extracted from the understanding 
of its vibration modes. The extraction of information from other solar-like stars has 
started to occur and the sustained modeling efforts by several groups has yielded many 
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asteroseismic studies (e.g. 


Thoul et al. 


2003; 


Eeeenbereer et al. 


2004 




Midio & Montalban 


2005 




Carrier et al. 


2005 




Eeeenbereer et al. 


2005 




2008 




Mosser et al. 


2008 


). This has 



positively confirmed the promise of asteroseismology: to dramatically improve knowledge 
about stellar parameters. This prospect, especially the promise to determine stellar radii 
within a 3 % error, is very attract ive for the exoplanet groups involved in the Kepler mission. 



For example, 



Baines et al. 



fl20o"l showed how the interferometric radius measurement of 



HD 189733 has dramatically improved and clarified the value for the radius a nd mean 
density of its planetary companion. We expect radius measurement from asteroseismology 



to have similar effect. An early indication o 
characterize exoplanet hosts was provided by 



the power of Kepler asteroseismol ogy to 



Christensen-Dalsgaard et al. 



mm . 



Photometric measurement of planetary transits, given that a precise modeling of the 
limb darkening is achieved (Prsa et al. 2010, in preparation), permits to identify four 
parameters of the star-planet systems: its orbital period Ilorb; its inclination i toward the 
observer and the fractional radii of the star and planet, which are defined as 

-Ra -Rb 

ta = , and = — , (1) 

a a 

respectively. We clearly see that the knowledge of the stellar host radius i?A and its 
associated error directly leads to a determination of the desired planet radius -Rb? as well 
as the orbital semi-major axis a. Further details on the influence of the precision of stellar 
parameters o n the inference of physical properties of planetary system were provided by 



SouthworthI (120091 ). One of the goals of the present study is to test the rehability of 
asteroseismology in determining the radius of stars by comparing the results with available 
interferometric measurements. 

This work presents an objective procedure to extract information from solar-like 
oscillations. The procedure, named SEEK, determines fundamental stellar parameters with 
the use of quantities inspired by the p-mode asymptotic asteroseismic relations, namely the 
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large and small separations, lS.v and (5z/, respectively (lTassoullll980l ). It can also make use 
of standard astrophysical input, the parallax vr, the V band magnitude and the reddening 
Eb-V: as well as input derived from atmosphere modeling, the effective temperature T^s, 
the gravity \ogg and finally the iron to hydrogen ratio [Fe/H]. Therefore, SEEK is not a 
purely asteroseismic procedure that only uses the frequency spectra of a pulsating star 
to constrain its fundamental parameters, but rather a hybrid procedure that uses both 
traditional and seismic input. 

The procedure is part of the Kepler Asteroseismic Science Consortium (KASC) pipeline 
that provides in priority stellar parameters to be used by exoplanet seekers of the Kepler 
team. Its aims are also to analyze every star in which Au and/or has been fixed and to 
provide an extended and homogeneous overview of the disk population observed by Kepler. 
A description of SEEK and of its main characteristics is made in the following sections. We 
tested SEEK on a set of well-studied stars. Of course the Sun is the golden standard, but we 
also tested the procedure against stars with diameters determined through interferometric 
measurements and/or with masses derived from Kepler's third law, thus testing our 
procedure against the closest we have to absolute quantities. The stars studied comprise 
/3 Hydri (HD 2151), r Ceti (HD 10700), Procyon A (HD 61421), r] Bootis (HD 121370), 
a Centauri A and B (HD 128620 and HD 128621), and 70 Ophiuchi A (HD 165341). 



2. THE SEEK PROCEDURE 

2.1. Stellar models 

The SEEK procedure makes use of a large grid of stellar models computed with the 
Aarhus Stellar Evolution Code (ASTEC). It compares an observed star with every model 
of the grid and makes a probabilistic assessment, with the help of Bayesian statistics, about 
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the structure of that star. 



Details of ASTEC and its input physics were described by IChristensen-Dalsgaard 

by A STEC, we used, for all our computations: 



( l2008al ). Out of the possible options offered 



the OPAL equation of st ate (ISwenson et al 



1996 



Alexander opacity tables (llglesias fc Rogerslll996 
to element ratios in the metallic mixture of 
treated with the mixing-length formulation of 



along with the OPAL plus Ferguson & 



Alexander fc 



erguson 



Grevesse &: Sauvall (119981) . and convection 



1994 ). the element 



Bohm-Vitensd (119581 ): the mixing length to 



pressure scale height ratio a, characterizing the convective efficacy, was treated as a variable 
parameter in the SEEK fits. Neither diffusion nor overshooting was included. 

A brief introduction to a beta version of SEEK was p resen ted along with a series 



Stello et al. 



( I2OO9I ). The beta version was 



of tests using synthetic observations in § 5 of 
built with an ad hoc probabilistic approach that is quite different from the one used in 
the present version. There is a major difference in the way errors on stellar parameters is 
obtained. This difference lies in the fact that the new grids have been extended and are 
dense enough to permit the use of Bayesian statistic. While SEEK yields similar results 
with both approaches, the Bayesian method used here is more rigorous, more robust and 
more comprehensive than the technique developed for the beta version of SEEK. 

The core of SEEK is a grid of models used to fit observations. We have calculated 
7,300 evolution tracks or 5,842,619 individual models. Each track calculation was 
stopped after the track had reached the giant branch or an age of r = 15 x 10^ yr. 
The tracks are separated into 100 subsets using a different combination of metallicity 
Z, initial hydrogen content Xq and mixing-length parameter a. These combinations are 
separated into two regularly spaced and interlaced subgrids. The first subgrid comprises 
tracks with Z = [0.005,0.01,0.015,0.02,0.025,0.03], Xq = [0.68,0.70,0.72,0.74], and 
a = [0.8,1.8,2.8] while the second subset has Z = [0.0075,0.0125,0.0175,0.0225,0.0275], 
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Xq = [0.69,0.71,0.73], a = [1.3,2.3]. Every subset is composed of 73 tracks spanning from 
0.6 to 3.0 Mq. The spacing in mass between the tracks is 0.02 Mq from 0.6 to 1.8 Mq and 
0.1 from 1.0 to 3.0 Mq. A shce of the grid is presented in Figure reffig:shce for Z = 0.015, 
Xq = 0.72 and a = 1.8 in a \ogg — logTefr diagram. 

The structure of the SEEK procedure is quite flexible. It will be easy in future to 
expand the grid with regularly or irregularly spaced tracks. In addition, we have started 
to test the possibility of interpolating models along the age-mass plane of the grid and 
between grid values of Z and Xq. The results are quite promising and it seems that the 
grid is dense enough with regards to these parameters to realize precise interpolation. 
However, interpolation tests in the a direction were not convincing and it seems that the 
grid would have to be refined significantly with respect to that parameter to permit precise 
interpolation. 

Some care has been required in defining the parameters of the grid, including the 
specific range and spacing we have chosen. Several factors influenced the definition of 
the compositi on. Recent analyses have led to substantial revisions of the solar surface 



composition ( lAsplund 



2005 



Grevesse et al. 



2007 



Asplund et al. 



20091 ): since the solar 



composition is used as reference to relate the observed [Fe/H] to Z/X this potentially 
affects the assumed stella r composition. Here we have chosen to use the older value of 



Grevesse fc Sauvall ( 119981 ). This can be justified based on the fact that solar models using 



the old values are closer to the helioseismic inferences than t 
abundances (see, for example 



Basu et al. 



120071 : 



lose u sing the more recent 



Basu &: Antia 



20081 ) . We also emphasize 



that for most stars, atmospheric parameters in the literature are still only available with 
older metallic mixtures. Homogeneity of our analysis requires the use of these mixtures. 
The metalhcity range used in the grid, —0.61 < [Fe/H] < 0.20, covers most (> 90%) of 
the stellar population of the disk if the Hipparcos color-magnitude diagram is taken as a 
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reference (lReidlll999l : IChabrierll200l[ ). This should permit the analysis of a large majority 
of the Kepler stars. 

The careful reader will also remark that relatively high values of Yq = 0.2713 and 
Zq = 0.0196 have been used in the standard definition of [Fe/H]. This value is used t o 



19981 ). 



calibrate solar models of ASTEC to the right luminosity ( IChristensen-Dalsgaard 
This value is used in all conversions done here. 

The mass step of 0.02 Mq from 0.6 to 1.8 Mq has first been chosen to be small enough 
to permit interpolation in the log g — log T^s diagram. However, we realized that this small 
step in mass has the virtue of making the grid sufficiently dense that we do not have to 
rely on interpolation to make our probabilistic estimate of the stellar parameters. The 
lower limit to the range results from the fact that molecules start to play an important role 
in these stars and that our evolution code does not include them. This should not be a 
handicap to SEEK since the low absolute magnitude of the stars not included in the grid, 
with M < 0.5 Mq, makes them improbable Kepler targets. The upper limit of 1.8 Mq is 
approximately the highest mass of the grid at which a model on the main sequence has 
at least 1% of its radius as a convective envelope. We wanted to keep a high density of 
models up to that limit where it is most likely to observe stars with stochasticall y excited 
oscillations. The models at higher masses are put as a safeguard permitting us to pin down 
a star with a very thin convective zone showing solar-like oscillations. We also remark 
that typically models with higher a and/or higher Xq and/or higher Z will tend to keep a 
convective envelope at slightly higher masses. Similarly, these models develop a convective 
core, leading to the typical hook path in the H-R diagram at higher masses. 

The grid mid- value of a = 1.8 was chosen in accordance with the solar calibration. 
The span of a is a trade-off between trying to get a wide range of values, from 0.8 to 2.8, 
and to limit the quantity of computed models in the grid. We had a bias towards testing 
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a wide range of convective efficiency since it is not clear how the mixing-length theory 
extends from the Sun to other main-sequence stars. However, this choice is not ideal since, 
as mentioned, the steps are too large to permit interpolation between the grid points. 

The role of SEEK is more to measure the size of the "valley of good solutions" produced 
by Kepler observations than it is to probe the deepest abyss of that valley. Specifically, the 
quality of a solution is measured by how well the computed observables in the grid, {qf}, 
match the observed values {qf^^}] this is quantified in terms of 

i=0 ^ « / 

where cxj is the estimated error for the observation q°^^ and n the number of observables. 
SEEK does not claim to find the best model for an observed star; its aim is to draw the 
contour of good solutions which is located around Xr = X^/P ~ 1, -P being the number of 
degrees of freedom of the problem. SEEK also outputs the stellar parameters with reliable 
error bars. We do not want to underestimate the errors by restricting the size of our grid. 
The approach is especially well suited to our problem since our solutions coming from the 
observation of Au are always heavily degenerate. This means that we do not expect a single 
and small region of the parameter space to yield a result that is much better than any other 
in terms of x^- Instead we have, as illustrated in Figure [2] for the star a Centauri A, a large 
region of the parameter space with ~ 1 . In Figure [2], we present for drawing convenience 
log Xr as a function of the mass M and the normalized age defined as 

= —T-. ^-777^ T , (3) 

mm(rRGB, 14Gyr) 

where r is the age of a model on a specific track and trgb is the age of the model 
following the same track when it reaches the red-giant branch (RGB). In Figure HI this 
parameterization is useful for visual purposes as it permits a smooth transition of the 
values from the lower-mass models not reaching the RGB to the models reaching it. The 
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exact position of trgb on the RGB is not critical since the evolution time scale on the 
RGB is much faster than on the main sequence. We stress that this parameterization was 
used for cosmetic purposes and is not involved in the SEEK computations. The models 
showing higher Xr to the left of the valley are caused by the appearance of convective cores 
in heavier models; this makes the automatic computation of the small separation somewhat 
difficult and creates these artifacts. 



Two works presenting pipe l ines t o do asteroseismo . 



been published, IMetcalfe et al. 



(120091) and 



Basu et al 



ogy o f solar-like stars have recently 



teOlOl ). In the former, the authors 



combine a pipeline using a genetic algorithm to find the best model in the parameter 
space and a local analysis, based on linearizatio n arou nd this solution, to refine the fit and 



estimate the error in that solution. 



Basu et al. 



( I2OIOI ) find the best model with a grid of 



stellar models, not unlike SEEK but using a smaller grid, while the errors are determined 
with a frequentist approach, using a Monte Carlo simulation of synthetic observations. 
We certainly recognize the desirability of a genetic algorithm, or a similar technique 
aiming at the absolute minimum in the parameter space, while also acknowledging the 
large computation power needed for suc h techniques. However, we note that a local error 



estimator used by 



Metcalfe et al. 



( I2OO9I ) assumes that the solution is linear within the 
error bars and that no nearby local minima can be deep enough to contribute to the error 
analysis; this can lead to underestimation of the errors. In our view, only an error estimator, 
Bayesian or frequentist, looking at various models around the best solutions may provide 
fully reliable results. This is particularly true in a problem, as the present, where the valley 
of solution is extended. This preference for Bayesian statistics has motivated our use of it 
in the development of SEEK. 
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2.2. The Observables 

The grid permits the mapping of the model physical input parameters p = 
{M, T, Z, Xq, a} into the grid of observable quantities = {Az/, 6^, Tgfj, log (7, [Fe/H], My, ...}, 
defining the transformation 

= /C(p). (4) 
We use these quantities and compare them to the actual observed quantities q°^^. 

Some of the observables used for the subsequent fit are easily extracted from the grid 
like Teff, logs', [F^/H]. However, My, Az/, and 5z/ need some extra attention. We compute 
absolute magnitude My , in the Johnson V b a nd, o f our model by applying the bolometric 
correction -Bc(^) of the 



VandenBerg &: Clenj (|2003[ ) tables to our luminosity L, 



where Mboi,© is the 



My = Mbo1,0 - 2.5 log — - B,{V), 



(5) 



VandenBerg fc ClemI ( l2003l ) prescription of 4.75 for the Sun's bolometric 



magnitude. We can then compare this value to its observed counterpart defined as 

M^^' = 5log7c - {V - A^) + 5. 

The conversion is done using the following relation between the absorption Ay and the 
interstellar reddening Eb-v'- Ay = 3.1Eb-v- 



We computed p-mode freq uencies for all models on the m ain sequence and subgiant 



(6) 



branch using the ADIPLS code (jChristensen-Dalsgaard 



2008bl ): the models on the red-giant 



branch are not included in this asteroseismic study. The large separation Au is computed 
with modes of / = only, 

Az/ = (z/„,o - i^n^i,o) , (7) 
while is the combination of modes with Z = 0, 2 only. 
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here Un^i is the frequency of the pulsation mode with angular degree / and radial order n. 
We know that these asymptotic values change with the radial order n in models. This is 
especially true in the case of Su. To compare observed separations with model values we 
have computed the average separations with up to 8 modes having different n and centered 
on a predetermined value of z/centrai- The quantity z/centrai is related to th e maximum powe r 



Stello et al. 



2009 



i/jnax seen in the Fourier transform of light curves of solar-like stars (see 
for comments on t'max)- Thus, in our grid the large and small separations are functions of 
the observed Uuiax- 

The computation of the average large and small separations is done by finding the 
maximum of the cross-correlation function of the frequencies selected around a suitable 
value t'centrai ncar z/max- Figure [3] shows how the small separation is computed for a specific 
model of the grid. In panel (a), the function ipQ corresponding to the / = eigenfrequencies 
is shown in blue and the function tp2 for / = 2 in red. These functions are convolved with 
the triangular distribution T of panel (b). The base of the triangular distribution is one 
fifth of the smallest distance between any combination of the / = eigenfrequencies. The 
convolution 



Di{x)= / ipi{s)T{x - s)ds 2 = 0,2 (9) 

J — oo 

is presented in panel (c) while the resulting cross-correlation, 

/oo 
Dois)D2ix + s)ds , (10) 
-oo 

is in panel (d). The maximum of the correlation is at Su. It is easy to see that the 
extra / = 2 mode around u ~ 750 /iHz is only producing a small bump in the correlation 
function and has no effect on the value of Su. For Az/, the location of the maximum of the 
auto-correlation [Dq * Dq), apart from the maximum at zero frequency shift, is taken. This 
technique is more robust than using the mean value of Equation ([7]) and especially ([H]) over 
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different values of n. The cross-correlation is not influenced by missing or extra eigenmodes, 
usually nonradial, produced by ADIPLS or by ambiguous determination of the radial order 

of a mode. The method is especially robust when it comes to computing Su for models 
having a convective core, or for models with avoided crossing. Robustness is very desirable 
since it is not possible to control the ADIPLS output for the 6 million models included in 
th e grid. 

In the present version of SEEK, we computed all adiabatic modes with / = 0, 2 in 
a fixed range of the dimensionless frequency, 10 < < 2800, where a is related to the 
angular mode frequency u by 

(") 

where G is the gravitational constant. This range makes it possible to cover observed 
solar-like pulsation in known stars on the main sequence and subgiant branch. These 
pulsations cluster around th e value of I'raax. which tends to follow the scaling relation of 



Kjeldsen fc Bedding 



(119951 ) 



M 

To accelerate the SEEK procedure and for models showing eigenmodes in these region, 
we precalculated 25 different sets of Au and for values of t'centrai equally spaced on 
a logarithmic scale, ranging from the Nyquist frequency of the short cadence of Kepler, 
8330 /iHz (corresponding to a period of 2 min) to 27.8 /iHz (corresponding to a period of 
10 h). (For any given model the actual range in t'centrai was restricted to values corresponding 
to acoustic modes for that model.) If, for example, we want to compare the grid value of 
5z/ with the solar value of Suq at z/max,o ~ 3333 /iHz, our grid automatically chooses the set 

'^'^( '^central) which is the closest to Z^max,© S't ^central = 3220 /iHz. 
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BAYESIAN APPROACH 



The Bayesian stat istics method of SEE K a nd the notation used he r e wer e inspired 
mainly by the work of 



J0rgensen &: LindegrenI ( 120051 ) . These 



Pont fc Eved (120041 ) and 
investigations compared the Bayesian approach to other means used for th e determ i nation 
of stellar age. We refer the reader to these papers and to the textbook by 



Gregory 



(120051) 



for technical questions involved in the use of Bayesian statistics. 



Let us first define the maximum likelihood function 



vi=0 



(13) 



with defined in Equation A maximum-likelihood estimate of stellar parameters can 
be obtained by finding the maximum of C, which, in case of Gaussian errors, is equivalent 
to minimizing x^- It can be argued that a maximum likelihood estimator is often enough 
to estimate stellar parameters. However, the highly non-linear mapping function K, makes 
it necessary to provide a set of priors, especially if a realistic estimate of the errors on the 
stellar parameters is to be made. The prior can be seen as a specific weight attached to 
every point of the grid determining the probability of tha t point to be observed. The ni ost 



obvious example of the utility of priors is well studied in iJeirgensen fc LindegrenI ([20051) 



It arises when one has to choose between different stellar models at different evolutionary 
stages but showing the same observables. In their example, one model is slowly evolving 
on the main sequence while the other is rushing up the subgiant branch. An experienced 
astronomer would naturally choose the slowly evolving star as the most probable model 
even if the resulting £ value is the same for both models. For the astronomer using Bayesian 
tools, a larger prior weight fo{p) given to the slowly evolving model automatically makes it 
a better candidate because the resulting probability, or posterior / defined at each point j 
as 

P oc gC^ , (14) 
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is bigger for that model. The grid overall posterior can be described as a density function, 

f{p) oc fo{p)C{JC{p)) . (15) 

iFiom the grid presented in the previous section, the computation of C at every point 
is straightforward. On the other hand, the definition of the priors / requires some insight 
in the problem. 

The most natural way to write the prior density is as a function of p, the model 
parameters, 

/o = vI/(r)$(Z|r)C(Xo|Z, tK{M\Xo, Z, r)/3{a\Xo, Z, r, M) (16) 

where ^1/(t) is the star formation rate (SFR) through time, and $(Z|r) and ({Xq\Z,t) 
describe the metallicity and initial hydrogen mixture as a function of age and as a function 
of age and metallicity, respectively. The initial mass function (IMF) C,{M\Xo, Z, r) is 
function of the element mixture as well as a function of time, and the prior related to the 
mixing-length parameter /3(a|Xo, Z, r, m) can also depend on the other stellar parameters. 

The distribution and the correlations can be built with the help of assumptions made 
from observations and from stellar and galactic models. In our case, since we are mainly 
interested in the Kepler satellite field of view which is well documented through the Kepler 
Input Catalog (KIC), we could use Equation @ to map a color-color distribution, or any 
other relation found in the KIC, into the prior of Equation (|16p. Ho wever, we shall keep here 



a more conservative approach, or as 



J0rgensen fc Lindegrenl (j2005[ ) put it, a non-committal 



approach and assume that r, Z, Xq, M and a are independent. Thus we write 

/o = vI/(r)<|.(Z)C(Xo)e(M)/3(«) , (17) 

and also take "^{r), ^{Z), C(-^o)5 and f3{a) to be flat. This means that the prior density is 
constant through these dimensions. We only make an assumption on the IMF and choose 
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the so-called IMFl model from IChabrierl (120011 ) where 



0.019 if M < l.OMo 

m) = < ~ (18) 

0.019 if M > l.OM©. 

In a grid, flat priors can be easily illustrated by equally spaced points, every point 
having an equal probability 1/A^. Of course, if the grid is not regularly spaced the weight 
assignment can get a bit more complicated. For a non-flat prior, we could distribute the N 
grid points so that they follow the density prescribed for example by Equation f[T5]) and 
give every point a weight of 1/A^ and meet the non-fiat IMF prior prescription. However, 
and this is especially true in the r direction, we have an irregularly spaced grid. Also, we do 
not wish to interpolate our grid for the sake of precision and for the same reason it is not 
necessarily reasonable to make the grid sparser in the region where the prior is relatively 
small. A star observed in that sparser region, even if less likely, would not be analyzed in 
the same detail by SEEK. 

We have solved this problem by dividing the parameter space into a number of small 
5-D bins of (penta-)volume containing each a small number of point > 1. Then, a 
prior weight /q, where 

r,= [ vI/(r)$(Z)C(Xo)e(M)/3(a)rfr , (19) 

is distributed through the point contained in the bin such that the weight /g^. given to 
model k is inversely proportional to the evolution speed of that model, 

/ofc = ^^^/o where k = 1,2,3, ...,n\ (20) 

Here Ar^ is the time taken by model k to evolve to the next point in its evolution track. 
This forces the slowly evolving models of the bin to have larger priors. These two equations 
are central pieces of SEEK. They ensure that the probability given to a model reflects the 
size of the parameter space where no other models are computed around it. These empty 
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regions of the parameter space will be seen as gaps in some of the probability distributions 
of the next section. 

4. RESULTS 

In practice, once / is fixed, the computation of the posterior probability Q on an 
arbitrary set of parameters is easily done. In one dimension, the grid is sliced along the 
desired parameter and all models entering a bin from s to s + As are added 

K 

g{s,) = Cj2fis.) , (21) 

k=0 

where the f^{si) are the K points of the grid lying from Sj to Sj+i = Sj + As, and C 
is a constant ensuring that the total probability of the problem is 1 (see Equation [22]) . 
The results of this computation is a probability distribution G{s) that is best seen as a 
histogram. The results for four parameters i?, M , r, and Teg, computed for 70 Ophiuchi A, 
are presented in Figure |H For convenience, there and in the following figures, we have 
normalized the histograms so that max(^(s)) = 1. 

Another important quantity is the running integral of Q seen as the solid red line in 
Figure |H The running integral over the entire parameter space is the total probability of 
the grid. This should be 1 since we assume that the observed star lies within the grid. We 
use that relation to fix the constant C in Equation (121 p 

Y,g{s,) = 1 . (22) 

i 

Similarly, we can integrate over a section of the distribution to get a 1, 2, or 3 a probability 
range for s. The error interval [s*^^\ s*^^^] output by SEEK is by default at the 1 a level and 
is determined based on the probability integral of Q. Specifically, 

/ g{s)ds = 0.1585 , / g (s) ds = 0.8A15 , (23) 

J So So 
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where sq is the initial value of the parameter, and the integrals are represented by suitable 
sums, as in Equation fl22]) . (Note that 0.8415 — 0.1585 = 0.6830 corresponds to 1 a for a 
normal distribution.) The solution SEEK outputs is the mid-point (s'-^^ + s*^^^)/2 of this 
interval. In this way the error bars are symmetric; this would not be the case if the median 



had been chosen. This 1 a section of the running integral is depicted in blue (between the 
dots) on the running integral of Figure HI 

These histograms along with the net numbers output by SEEK are useful to 
interpret the results. In the case of 70 Oph A, the solutions R = 0.86 ± 0.01 [Rq] and 
Teff = 5306 ± 36 [K] are well constrained and quasi Gaussian. The histograms visually 
confirm the validity of the numerical solutions. The result for the mass M = 0.90 ± 0.04 Mq 
reveals some of the grid's limitation as there are gaps in the otherwise well defined structure 
of Q{M). These gaps are artifacts of the 0.02 Mq mesh used in the grid below 1.8 Mq. 
However, since the prior density is properly weighted through Equations (1211) and (1221) . the 
overall value of the mass, with its error bars, is correct. For the age, r = 8.65 ± 3.64 Gyr, a 
large error is revealed by the output and the histogram shows the probability distribution to 
be populated all along r from 2 Gyr all the way to 13 Gyr. In that case a more reasonable 
answer is a lower limit of r > 2 Gyr. The histogram for r also shows gaps that are due to 
the finite resolution of the grid. 

Table [T] shows our selection of solar-like oscillators and observables found in the 
literature. These observables are processed by SEEK and produce the outputs in Table [2] 
and Table [31 In Table [21 we compared SEEK results with direct measurements of radius via 
interferometry and of mass via Kepler's third law. For SEEK, the fits have been made for 



value s*^™'^ with 




(24) 
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two different posterior sets, the full set 

q^'' = {Au, Teff, log^, [Fe/H], V, vr, Eb-v} , (25) 

and one using only a selected subset of the same observables, 

gg'^ = {Az/, Teff, [Fe/H]}. (26) 

The subset may be more representative of what we can expect from a typical Kepler 
observation when no ground-based follow-up has been done for a star, and where the 
asteroseismic data have not yielded a reliable determination of 5v. This can occur in a 
stochastically oscillating star when the width of observed modes is larger than the small 
separation itself. In that case, we have Az/ from the satellite observations, while Teg and 
[Fe/H] can be obtained from the KIC. 

The results of Table |2] show that the 3% promise on the radius precision mentioned in 
the introduction is generally reached with both the full set and the subset of observables. 
More striking and certainly more important than the precision itself is the very good 
accuracy reached by SEEK. It is on this basis that our method, or any other technique, 
should be judged, provided that independent measurements are available. SEEK has pinned 
down all stellar radii at the 3% level when the extended set of observable was used. The 
same level of accuracy is generally reached if the subset of observables (Eq. [2^ is used. 

We note that in well-posed problems, where parameters are actually strongly 
constrained, the value of C is much larger than /o close to the best solution. It is only 
when we get away from the best parameter values that the prior gets more important. The 
details of the priors used here only have an influence on the wing of the Q distribution. 
Thus, the choice of priors has an effect on the errors, especially if the errors are large, but 
little effect on the solution itself. 
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DISCUSSION 



We clearly need to understand the exceptions to the general success of SEEK; in 
particular, we note that /3 Hyi was a bit further from the direct observations with an offset 
of 6%. An extra tool offered by SEEK can be used to study the P Hyi radius offset when 
compared to direct measurements. Figure [5] shows a two-dimensional projection Q{R, Z/X) 
of the posterior computed for /3 Hyi with only Au and T^s as input observables. The 
metallicity is not included as an observable in the computation as a trick to have a clearer 
view of the correlation between that observable and R. Several possibilities are offered 
to explain the discrepancy between SEEK's radius R/Rq = 1.92 ± 0.05 and the direct 
measurement of 1.814 ± 0.017. The first possibility is that the observed Z/X is too high. 
Reducing [Fe/H] by ~ 0.2 dex, to Z/X ~ 0.0163, would reconcile the asteroseismic and 
interferometric radii. This can be concluded from a visual examination of Figure [51 In fact, 
lowering the metallic ity by this amount would put it back to the level measured in /3 Hyi by 



Dravins et al 



( 119981 ). We have made an a post eriori test using the value of metallicity and 



temperature published by 



Dravins et al 



instead of the more recent values of Table [T] (Ida Silva et al.ll2006l ). These inputs processed 



(Il998h. T.(T = 5800 ± 100 K. and [Fe/H]= 0.2 ± 0.1 



by SEEK place the radius RlR(^ = 1.85 



observations. We note that 



Dravins et al. 



0.03 in much better agreement (2%) w' i th the 



(1l998h used the older iNoels fc Grevessel (11993[ l 



metal mixtures in the computation. This underlines that accuracy can only be reached if 
the input metallicity is selected with some care. 

The offset in radius for P Hyi can be explained otherwise. The extension of the 
grid to lower values of Xq would to some extent reconcile SEEK and the interferometric 
measurement. It can be seen in Figure [5] that the left-hand side of the correlation function 
Q is in fact the edge of the grid where Xq = 0.68. The extension of the grid would likely 
expand the existing ridge of solutions to the left in Figure [5l at lower radii. A simple 
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eye-ball estimate in the figure can be made to indicate that the combination, Z ~ 0.0165, 
Xq ~ 0.64 would yield the right metallicity and the right radius. However, it seems that 
a hydrogen content that low is not reasonable and perhaps this possibility should not be 
considered seriously. The effect of the hydrogen content, or should we say helium content, 
on the radius is in fact the main factor compromising the precision of R in SEEK or any 
other technique using the large sepa ration to determine the radii. The error on the radius 
is influenced by the size of the grid in the Xq direction. The imperative we had not to 
underestimate the errors on stellar parameters influenced our choice of having a large range 
of values for that parameter. It also means that raising the precision on the metallicity 
ratio ZjX to fix the radius with better precision has some intrinsic limitations, unless the 
helium abundance can be fixed independently. We also note that the gaps between the 
islands of Figure like the one for the ID histograms, are caused by the finite resolution of 
the grid. The gaps in this figure also reflect that if the grid is regularly spaced in Xq and Z 
it is not m. ZjX since the gaps are not regularly spaced along that direction. 

One other result obtained by SEEK is worth examining in more detail. For Procyon A, 
if the subset of observables is used to do the fit, the values of M and R are less precise, 
which is in accordance with expectation, but more accurate than the case where all 
observables are used. This could reveal some inconsistency in the observables as well as 
some of the limitations of the SEEK procedure. Indeed forcing SEEK to look for a model 
of Procyon A with the right astrometry, V = 0.363 and vr = 134.07 mas, reduces the 
accuracy of the results, compared with using only Tes, [Fe/H] and Au. This discrepancy 
could come from inhomogeneity between the tool used in the spectroscopic study giving 
the observables and in our models . The metallicity and tern perature taken as input for 



Procy on A are from 3D mo 



of the 



d eling (lAUende Prieto et al. 



20021 ) . while we used the conversion 



VandenBerg fc ClemI ( 120031 ) tables using ID models to get our value of My. Using 



older atmospheric parameters, or a f ully consistent conversion, could perhaps improve the 
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fit. 

We also remark that the behavior of our results for the mass of a Cen B presents 
similarity to those of Procyon A. The results are slightly more accurate if V ^ ir and Su are 
dropped from the fit. It seems that this reveals some of the limitation of our grid when 
extremely precise measurements are available for a star. The small box defined within 
~ 1 — o" of all a Cen B observables only includes a few models. This number is too small and 
can only produce a weak probabilistic assessment over the value of the mass. Computing 
more models to cover with more precision the parameters space of a Cen B would most 
likely resolve the small deviation over the value of the mass. This underscores the fact that 
modeling dedicated to individual stars becomes useful once classical observables of a star 
are very precise. 

In any case, Table |2] shows that SEEK is certainly capable of pinning down the radius 
of solar-like oscillators with a minimum effort and that it is also able to fix the mass of these 
stars with an accuracy not attainable when only non asteroseismic inputs are available. 
This improvement in accuracy is especially potent in case the distance to the star is not 
known. 

In Table El all primary stellar parameters p of the modeled stars are shown as well as 
some selected secondary values. We address here the question of the age of stars fixed by 
SEEK, and especially the age of the Sun, since it is the only star where an independent 
determination of the age is available. A fairly good accuracy, 3.9 6=b0.41 Gyr (13%), 
compared to the meteoritic value of 4.57ib0.02 Gyr (Wasserburg, in 



Bahcall &: Pinsonneault 



19951 ). is reached in the case where Au and are known. A look at the histogram in the 
top panel of Figure |6] confirms visually that the solution is well constrained. This is quite 
promising in the case where the small separation is known. It seems that we can expect a 
precision of 5 to 20% for the age when Az/, Tes, [Fe/H], and especially Su are known with 
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reasonable precision. 

If only the three observables in the subset are used, the result varies from star to star. 
For the Sun, the lower panel of Figure [6] shows for the age two distinct and very probable 
solutions. This specific example shows how useful the histograms are since the output of 
Table [3] does not reveal the presence of two solutions. One thing not shown here, but that 
we were able to see from the computation of the probability distribution ^(r, a), is that 
the two islands of solution for the age are at two different value of a, the young solution 
corresponding to the standard value of a = 1.8 and the older one to a = 2.8. As seen in the 
lower panel of Figure M higher values of a are acceptable if only the subset of observables is 
used. 

Figure [7] illustrates another type of effect on the age when the subset of observables 
is used. The precision on the age of /3 Hyi is degraded by 20%, but the overall answer, or 
accuracy, is the same. As mentioned previously in the description of Figure HI 70 Oph A 
shows a third type of scenario where not much can be concluded on the age of the star as 
the posterior is showing solutions at almost any age. Since the relation between the age 
and the observables is highly non-linear, it is very difficult to estimate beforehand what will 
be the precision obtained on the age depending on the type of star and/or on the precision 
of the observables, especially when 5z/ is not known. 

We also note that any age determination based on stellar modeling is susceptible to 
changes in basic aspects of the modeling. In the solar case helioseismic fits to low-degree 
modes result in a solar age very close to the meteoritic value, for models including diffusion 
and settling of heli um and heavy e l emen ts, and when the 'old' solar composition is used. 



whereas use of the 



Asplund et al. 



Christensen-Dalsgaard 



2009 



(120091 ) composition induces a significant shift (see 



for a review). Also, the age determination based on fits to the 



Kepler data for HAT-P-7, which has a convective core, was understandably rather sensitive 
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to th e extent to which convective-core overshoot was included (jChristensen-Dalsgaard et aL 



20101 ). Such potential systematic effects must clearly be taken into account in the 
interpretation of results of fits such as those carried out with SEEK. Of course the long-term 
goal is to reduce these effects through an improved understanding of stellar structure and 
evolution based on more detailed asteroseismic investigations. 



6. CONCLUSIONS 

We have shown the details involved in the SEEK procedure and we have tested the 
validity of the approach with all independent measurement known to us. This is the first 
step before Kepler ob servations of solar-like stars, now becoming available in large quantity 



( IGilliland et al. 



20101 ). are analyzed. We hope that SEEK will be able to process these 
observation with ease and lead to the publication of a catalog that includes a homogeneous 
sample of stars revealing their radius, mass, age, and other stellar parameters. 

We expect that the results of SEEK can be a useful starting point in more detailed 
asteroseismic analyses in the, likely frequent, cases where extensive sets of oscillation 
frequencies are available. Such detailed investigations can further refine the basic stellar 
parameters, including the age which, as shown above, is sometimes not significantly 
constrained by SEEK. Also, detailed analyses of the frequencies are likely to uncover 
evidence for the need for improvements in stellar modeling. 

The SEEK model grid extends to the base of the red-giant branch and hence SEEK can 
be used to determine properties of stars on both the main sequence and the subgiant branch. 
In fact, the results obtained for the subgiants 1] Boo and /3 Hyi indicate the potential 
of SEEK for the investigation of subgiants, common amongst the Kepler asteroseismic 
targets. On the other hand we acknowledge the substantially different diagnostic potential 
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between centrally hydrogen burning stars and subgiants. In particular, provides a direct 
measure of stellar age in the former case, while the diagnostic potential of 5v for subgiants 



is more subtle. In the latter case, detailed analysis o. 
may provide much more stringent constraints (e.g. 



individual frequenc ies of mixed modes 



Metcalfe et al. 



20 id ). These issues of 



optimizing the asteroseismic diagnostics certainly require further investigations. 

In parallel, further tests and development of SEEK are required. A first step will be 
to test the level of systematic errors that are introduced by the neglect of relevant physical 
effects in the stellar modeling. An obvious example is diffusion and settling which have been 
demonstrated to have a substantial effect on solar modeling, as tested with helioseismology. 
Also, effects of convective core overshoot could influence the results for stars slightly more 
massive than the Sun. New grids of models including such effects can then be computed for 
inclusion in SEEK. 

A related issue, with less obvious solutions, concerns the systematic effects of the 
frequency errors introduced by the failure to model properly the outermost layers of the 
star, which are known to dominate the difference between observed and modeled frequencies 



in the solar case (e.g.. 



Christensen-Dalsgaard et al. 



19961 ) and which also have potentially 



significant effects on the large and small separations. It is important at least to obtain 
an estimate of the extent to which these effects may influence the results obtained with 
SEEK. One may hope that improvements in the modeling of the relevant effects, combined 
with detailed analyses of Kepler data for a broad range of stars, will eventually allow us to 
reduce or eliminate such potential systematic errors. 
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Fig. 1. — Slice of the SEEK grid projected in the log (7 — logTcfj diagram for Z = 0.015 
Xo = 0.72 and a = 1.8. Three isochrones are also drawn. The computation was stopped 
after the star had reached 15 billion years or an arbitrary point on the RGB branch. The 
low-mass models are separated by steps of 0.02 Mq, the more massive ones by 0.10 Mq. 
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Fig. 2. — Slice of the solution valley obtained by SEEK for a Cen A at Z = 0.025, Xq = 0.70, 
and a = 1.8. Note the logarithmic color scale showing the value of Xr ^ function of M 
and tn (see Eq. [3l). The solution is centered on the constant radius line of R = 1.25 Rq, 
shown as a solid red line. The solid black lines are also radius isocurves. Models older than 
13 Gyr, on the upper right region of the figure, and the models on the ZAMS, with = 0, 
were excluded from the fit and given values of logXr — 6-0. The slight deviation in the 
constant-radius isocurves with R > 1.5 Rq and tn ^ 0.7 marks the position of the hook in 
the evolution tracks. 
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Fig. 3. — Description of the 6v computation. The functions i/jo,iIj2 corresponding to the 
computed frequencies for / = 0,2, shown in blue and red respectively in panel (a), are 
convolved with the distribution T in panel (b) to create the distribution Dq,2 in panel (c). 
The resulting cross-correlation function {Dq * D2) fixing the at its maximum is shown in 
panel (d). The extra bump of the cross-correlation is caused by the extra / = 2 mode at 
u ~ 760 /iHz and has no influence on the computed separation. 
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Fig. 4. — SEEK's posterior for 70 Oph A. All the observables of Table [T] are used in the fit. 
The probability distribution G{p) for R, M, t and T^s, (normalized to a maximum of 1, black 
histogram) is plotted along with the probability integral (solid red line). The delimitation 
of the 1(7 probability of the distribution, determined by the interval [s*-^-*, s^^-*] (cf. Eq. [2^ . 
is plotted between the blue dots. 
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Fig. 5. — Correlation between ZjX and R for /3 Hyi. The two thick sohd curves through 
peaks in the correlation function Q{R, Z / X) follow the constant metallicity ridges for Z = 
0.025 and 0.015. Along these curves the hydrogen mass ratio varies from the upper left peak 
to the lower right one as Xq = 0.68, 0.70, 0.72, and 0.74. A clear general trend is that R 
increases with Z and that it also increases with increa sing Xp. The solid h orizontal line is 



20061), with its error 



the 13 Hyi metalhcity ratio of Z/X = 0.0258 ± 0.003 (Ida Silva et al. 
bars as dashed lines, while the solid and dashed vertical lines show the direct interferometric 
measurement of R = 1.814 ±0.107. Note that the merging of the constant metallicity ridges 
at lower values is an artifact of the drawing procedures. 
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Fig. 6. — SEEK's posterior for the Sun. The top 9 panels represent the probabihty dis- 
tribution (normahzed to 1) as well as the probability integral for the model using all the 
observables of Table [H The lower panels are similar but use a subset of observables (Az/, 
Teff, and [Fe/H]). 
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Fig. 7. — Same as Figure [6] for [3 Hyi. 



Table 1. Observables 



star At' [mHz] Siy [fJ-Hz] i^max [mHz] T(,ff [K] [Fc/H] log g [cgs] V tt [mas] E^^y Reference 



Hyi 


57.2 ± 0.5 


5.3 ± 0.5 


1000 


5964 ± 100 


-0.03 ± 0.05 




2.8 ± 0.01 


134.07 ± 0.11 


0.01 


1 


a 


T Get 


169.6 ± 0.5 


12.7 ± 1.2 


4500 


5264 ± 100 


-0.5 ± 0.03 


4.36 ± 0.15 


3.39 ± 0.01 







2 


b 


Procyon A 


55.5 ± 0.5 




1000 


6514 ± 27 


-0.05 ± 0.05 




0.363 ± 0.003 


285.93 ± 0.88 





3 


c 


rf Boo 


39.9 ± 0.9 


4.0 ± 1.0 


750 


6007 ± 255 


0.09 ± 0.01 




2.68 ± 0.03 


88.17 ± 0.75 





4 


d 


a Cen A 


105.5 ± 0.5 


5.6 ± 0.7 


2410 


5847 ± 27 


0.24 ± 0.03 


4.34 ± 0.12 


-0.01 ± 0.01 


747.1 ± 1.2 





5 


e 


a Can B 


161.5 ± 0.5 


10.1 ± 0.6 


4100 


5316 ± 28 


0.25 ± 0.04 


4.44 ± 0.15 


1.35 ± 0.01 


747.1 ± 1.2 





6 


e 


70 Oph A 


161.7 ± 0.8 




4500 


5300 ± 50 


0.04 ± 0.05 




4.19 ± 0.014 


194.2 ± 1.2 





7 


f 


Sun 


134.8 ± 0.5 


9.1 ± 0.2 


3034 


5778 ± 20 


± 0.01 


4.44 ± 0.01 








8 


S 



References. — As tero seismology (1) Beddinge^ 



iKjeldsen et al, 



I I2OO5I : (7)[ 



Carrier &: Eggenberge: 



ai1l20o4 



(8) 



(2) Teixeira^^al 



Thiery et al 



I bood ; (3) lEggenberger et ahlbooj ; 



4lCarrier et 



r5 1Bouchv &: Carrie; 



jl bood : 



(6) 



Re fere nces. — Other observab les (a) da Silva et al 
\200i ({) lEggenberger et al.l[20o3 ; (g) \g 



revesse &: Sauval 199 



200^: (b ) Isoubiran et alj[l993 ; fcl lAllende Prieto et aljbooj : f d) llVIorossi et al.l[20oj : fel lporto de Mello et all 



3 



Table 2. Comparison With Direct Measurements 



R [Rq\ 



Star 


Selected 


All 


Direct /Interfcromctry 


Reference 


Selected 


All 


Direct/Kepler 's Law 


Reference 


/3 Hyi 


1.92 ± 0.05 


1.87 ± 0.01 


1.814 ± 0.017 


1 


1.27 ± 0.11 


1.16 ± 0.02 






T Cot 


0.79 ± 0.01 


0.81 ± 0.01 


0.790 ± 0.005 


2 


0.77 ± 0.03 


0.83 ± 0.01 






Procyon A 


2.04 ± 0.05 


2.12 ± 0.02 


2.048 ± 0.025 


3 


1.45 ± 0.10 


1.59 ± 0.05 


1.497 ± 0.037 


4 


rj Boo 


2.65 ± 0.08 


2.69 ± 0.06 


2.672 ± 0.028 


5 


1.62 ± 0.13 


1.71 ± 0.08 






a Ccn A 


1.23 ± 0.04 


1.25 ± 0.01 


1.224 ± 0.003 


6 


1.09 ± 0.09 


1.11 ± 0.01 


1.105 ± 0.007 


6 


Q, Ccn B 


0.87 ± 0.01 


0.87 ± 0.01 


0.863 ± 0.005 


6 


0.92 ± 0.04 


0.89 ± 0.01 


0.935 ± 0.006 


6 


70 Oph A 


0.86 ± 0.02 


0.86 ± 0.01 






0.89 ± 0.06 


0.89 ± 0.05 


0.890 ± 0.020 


7 


Sun 


1.02 ± 0.03 


1.01 ± 0.01 


1 




1.03 ± 0.08 


1.01 ± 0.02 


1 
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Table 3. SEEK Output 



Stellar Parameters 

Primary Secondary- 



Star Observables M [Mq] t [Gyr] Z Xq a. R [Rq] \o%{L/Lq) log g [cgs] Tgff [K] Y 



P Hyi 


All 


1, 


,16 ± 0, 


,02 


5.41 


± 


0, 


,89 


0, 


,018 


± 


0, 


,001 


0, 


,708 


± 


0, 


,022 


1 


,74 ± 


0, 


,60 


1 


,87 ± 


0, 


,01 


0.58 ±0.02 


3.964 


±0, 


,004 


5845 ± 38 


0, 


,274 


± 


0, 


,023 




Selected 


1. 


,27 ± 0, 


,11 


3.94 


± 


1 


,45 


0, 


,019 


± 


0, 


,002 


0, 


,710 


± 


0, 


,022 


1 


,69 ± 


0, 


,62 


1 


,92 ± 


0, 


,05 


0.62 it 0.05 


3.97 


±0, 


,02 


5959 ± 111 


0, 


,271 


± 


0, 


,024 


T Get 


All 


0, 


,83 ± 0, 


01 


0,75 


± 





,03 





007 


± 





,001 





730 


± 





,003 


1 


,30 ± 


0, 


17 


0, 


,81 ± 





,01 


-0.32 ± 0,02 


4.547 


± 0, 


OOfi 


5450 ± 34 





263 


± 





,004 




Selected 


0, 


,77 ± 0, 


,03 


6.31 


± 


2 


,18 


0, 


,007 


± 


0, 


,001 


0, 


,729 


dz 


0, 


,004 


1 


,30 ± 


0, 


17 


0, 


,79 ± 


0, 


,01 


-0.36 ± 0.03 


4.53 


± 0, 


,01 


5286 ± 71 


0, 


,264 


± 


0, 


,005 


Procyon A 


All 


1, 


,69 ±0, 


,06 


1.83 


± 


0, 


,24 


0, 


,021 


± 


0, 


,003 


0, 


,731 


± 


0, 


,012 


2, 


,42 ± 


0, 


,46 


2, 


,12 ± 


0, 


,02 


0.88 ±0.02 


3.991 


±0, 


,010 


6548 ± 37 


0, 


,248 


± 


0, 


,015 




Selected 




,45 i 0. 


, 10 




_(_ 





,64 


0. 


,018 


_|_ 





,002 


0. 


,710 







,023 




,93 it 


0, 


62 




,04 it 





,05 


O.S.i It U.Oii 


3.98 


_i_ n 

± U. 


,02 


DOlO It D4 


0, 




_|_ 





,025 


T7 Boo 


All 


1, 


,71 ± 0, 


,08 


1.90 


± 


0, 


,38 


0, 


,025 


± 


0, 


,001 


0, 


,729 


± 


0, 


,013 


1 


,32 ± 


0, 


,60 


2, 


,69 ± 


0, 


,06 


0.97 ± 0.02 


3.81 


± 0, 


,02 


6166 ± 80 


0, 


,246 


± 


0, 


,014 




Selected 


1. 


,62 ± 0, 


,13 


2.24 


± 


0, 


,62 


0, 


,025 


± 


0, 


,001 


0, 


,728 


± 


0, 


,013 


1 


,33 ± 


0, 


,64 


2, 


,65 it 


0, 


,08 


0.88 ± 0.08 


3.80 


±0, 


,02 


5974 ± 280 


0, 


,247 


± 


0, 


,014 ^ 

c 


a Cen A 


All 


1, 


,11 ± 0, 


,01 


4.95 


± 


0, 


,04 


0, 


,030 


± 


0, 


,001 


0, 


,680 


± 


0, 


,003 


1 


,80 ± 


0, 


17 


1 


,25 ± 


0, 


,01 


0.18 ± 0.02 


4.299 


± 0, 


,001 


5850 ± 34 


0, 


,290 


± 


0, 


,004 




Selected 


1, 


,09 ± 0, 


,09 


8.04 


± 


4, 


,05 


0, 


,030 


± 


0, 


,001 


0, 


,698 


± 


0, 


,020 


2, 


,25 ± 


0, 


,59 


1 


,23 it 


0, 


,04 


0.20 ± 0.03 


4.30 


±0, 


,01 


5850 ± 37 


0, 


,272 


± 


0, 


,021 


a Cen B 


All 


0, 


,89 ± 0, 


,01 


6.96 


± 


0, 


,03 


0, 


,030 


± 


0, 


,001 


0, 


,680 


± 


0, 


,003 


1 


,80 ± 


0, 


,17 


0, 


,87 it 


0, 


,01 


-0.32 ± 0.02 


4.621 


±0, 


,001 


5250 ± 34 


0, 


,290 


± 


0, 


,004 




Selected 


0, 


,92 ± 0, 


,04 


6.99 


± 


4, 


,45 


0, 


,029 


± 


0, 


,001 


0, 


,709 


± 


0, 


,030 


2, 


,23 ± 


0, 


,52 


0, 


,87 it 


0, 


,01 


-0.26 ± 0.03 


4.63 


±0, 


,01 


5309 ± 66 


0, 


,262 


± 


0, 


,020 


70 Oph A 


All 


0, 


,89 ± 0, 


,06 


8.66 


± 


3, 


,70 


0, 


,023 


± 


0, 


,003 


0, 


,716 


± 


0, 


,018 


2, 


,06 ± 


0, 


,34 


0, 


,86 it 


0, 


,01 


-0.27 ± 0.02 


4.626 


±0, 


,006 


5311 ± 65 


0, 


,261 


± 


0, 


,021 




Selected 


0, 


,89 ± 0, 


,06 


7.64 


± 


4, 


,54 


0, 


,022 


± 


0, 


,003 


0, 


,710 


± 


0, 


,026 


1 


,99 ± 


0, 


33 


0, 


,86 it 


0, 


,02 


-0.28 ± 0.03 


4.62 


±0, 


,01 


5304 ± 71 


0, 


,268 


± 


0, 


,029 


Sun 


All 


1. 


,01 ± 0, 


,02 


3.96 


± 


0, 


,41 


0, 


,020 


± 


0, 


,001 


0, 


,707 


± 


0, 


,011 


1 


,80 ± 


0, 


,17 


1 


,01 ± 


0, 


,01 


0.02 ± 0.02 


4.432 


±0, 


,006 


5753 ± 36 


0, 


,273 


± 


0, 


,012 




Selected 


1, 


,03 ± 0, 


,08 


6.80 


± 


6, 


,43 


0, 


,020 


± 


0, 


,001 


0, 


,729 


± 


0, 


,013 


2, 


,18 ± 


0, 


,53 


1 


,02 it 


0, 


,03 


0.03 ± 0.03 


4.44 


±0, 


,01 


5757 ± 39 


0, 


,261 


± 


0, 


,014 



